Search for: All records

Over the last decade, substantial progress has been made in understanding the topology of quasi-two-dimensional (2-D) non-equilibrium fluid flows driven by ATP-powered microtubules and microorganisms. By contrast, the topology of three-dimensional (3-D) active fluid flows still poses interesting open questions. Here, we study the topology of a spherically confined active flow using 3-D direct numerical simulations of generalized Navier–Stokes (GNS) equations at the scale of typical microfluidic experiments. Consistent with earlier results for unbounded periodic domains, our simulations confirm the formation of Beltrami-like bulk flows with spontaneously broken chiral symmetry in this model. Furthermore, by leveraging fast methods to compute linking numbers, we explicitly connect this chiral symmetry breaking to the entanglement statistics of vortex lines. We observe that the mean of linking number distribution converges to the global helicity, consistent with the asymptotic result by Arnold [InVladimir I. Arnold – Collected Works(ed. A.B. Givental, B.A. Khesin, A.N. Varchenko, V.A. Vassiliev & O.Y. Viro), pp. 357–375. Springer]. Additionally, we characterize the rate of convergence of this measure with respect to the number and length of observed vortex lines, and examine higher moments of the distribution. We find that the full distribution is well described by a k-Gamma distribution, in agreement with an entropic argument. Beyond active suspensions, the tools for the topological characterization of 3-D vector fields developed here are applicable to any solenoidal field whose curl is tangent to or cancels at the boundaries of a simply connected domain. 

Abstract. It is generally accepted that modeling Farley–Buneman instabilities requires resolving ion Landau damping to reproduce experimentally observed features. Particle-in-cell (PIC) simulations have been able to reproduce most of these but at a computational cost that severely affects their scalability. This limitation hinders the study of non-local phenomena that require three dimensions or coupling with larger-scale processes. We argue that a form of the five-moment fluid system can recreate several qualitative aspects of Farley–Buneman dynamics such as density and phase speed saturation, wave turning, and heating. Unexpectedly, these features are still reproduced even without using artificial viscosity to capture Landau damping. Comparing the proposed fluid models and a PIC implementation shows good qualitative agreement. 

We describe CATKE, a parameterization for fluxes associated with small‐scale or “microscale” ocean turbulent mixing on scales between 1 and 100 m. CATKE uses a downgradient formulation that depends on a prognostic turbulent kinetic energy (TKE) variable and a diagnostic mixing length scale that includes a dynamic convective adjustment (CA) component. With its dynamic convective mixing length, CATKE predicts not just the depth spanned by convective plumes but also the characteristic convective mixing timescale, an important aspect of turbulent convection not captured by simpler static CA schemes. As a result, CATKE can describe the competition between convection and other processes such as shear‐driven mixing and baroclinic restratification. To calibrate CATKE, we use Ensemble Kalman Inversion to minimize the error between 21 large eddy simulations (LESs) and predictions of the LES data by CATKE‐parameterized single column simulations at three different vertical resolutions. We find that CATKE makes accurate predictions of both idealized and realistic LES compared to microscale turbulence parameterizations commonly used in climate models. 

We numerically and theoretically investigate the Boussinesq Eady model, where a rapidly rotating density-stratified layer of fluid is subject to a meridional temperature gradient in thermal wind balance with a uniform vertically sheared zonal flow. Through a suite of numerical simulations, we show that the transport properties of the resulting turbulent flow are governed by quasigeostrophic (QG) dynamics in the rapidly rotating strongly stratified regime. The ‘vortex gas’ scaling predictions put forward in the context of the two-layer QG model carry over to this fully three-dimensional system: the functional dependence of the meridional flux on the control parameters is the same, the two adjustable parameters entering the theory taking slightly different values. In line with the QG prediction, the meridional heat flux is depth-independent. The vertical heat flux is such that turbulence transports buoyancy along isopycnals, except in narrow layers near the top and bottom boundaries, the thickness of which decreases as the diffusivities go to zero. The emergent (re)stratification is set by a simple balance between the vertical heat flux and diffusion along the vertical direction. Overall, this study demonstrates how the vortex-gas scaling theory can be adapted to quantitatively predict the magnitude and vertical structure of the meridional and vertical heat fluxes, and of the emergent stratification, without additional fitting parameters.